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Abstract 

We consider the effect of replacing in stochastic differential equations 
leading to the dynamical collapse of the statevector, white noise stochastic 
processes with non white ones. We prove that such a modification can 
be consistently performed without altering the most interesting features 
of the previous models. One of the reasons to discuss this matter derives 
from the desire of being allowed to deal with physical stochastic fields, such 
as the gravitational one, which cannot give rise to white noises. From our 
point of view the most relevant motivation for the approach we propose 
here derives from the fact that in relativistic models the occurrence of 
white noises is the main responsible for the appearance of untractable 
divergences. Therefore, one can hope that resorting to non white noises 
one can overcome such a difficulty. We investigate stochastic equations 
with non white noises, we discuss their reduction properties and their 
physical implications. Our analysis has a precise interest not only for 
the above mentioned subject but also for the general study of dissipative 
systems and decoherence. 



1 Introduction 

The aim of dynamical reduction models ^ — Q is to combine the Schrddinger 
evolution and the wavepacket reduction postulate into one universal dynamical 
equation, which is assumed to govern all physical processes. In this way, such a 
dynamics accounts both for the quantum properties of microscopic systems and 
for the classical properties of macroscopic ones. 

*e-mail: bassi@ictp.trieste.it 
^e-mail: ghirardi@ts.infn.it 
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This goal is achieved by adding to the Schrodinger equation new stochastic 
terms which induce a diffusion process — guided by a set of Gaussian white 
noises — of the statevector in Hilbert space: it is precisely this sort of "random 
walk" which is responsible for the localization mechanism. We will review these 
attempts in section 2. 

The main aim of this paper is to generalize the basic equations of dynamical 
reduction models to the case of general Gaussian noises, not necessarily white 
in time. There are three reasons for considering models of this kind: first of 
all, it is interesting to analyze if, and to which extent, the nice features of the 
reduction mechanisms proposed up to now depend on the white-noise character 
of the stochastic processes. Secondly, it has been argued several times that the 
stochastic processes should be related to physical fields — the most promising 
being the gravitational field; in such a case, the stochastic processes cannot 
be white, since white noises are never realized in nature. Finally, the third 
motivation to look for generalizations of the previous models has to do with the 
important problem of working out relativistic models of wavepacket reduction 
[|. In fact, even though it has been proved that models of this type share all the 
nice features of the non relativistic ones, they suffer from the serious drawback 
of inducing an infinite increase of the energy of physical systems. This infinite 
increase is basically caused by the local coupling between the quantum fields 
and the white-noise stochastic fields appearing in such theories. Accordingly, a 
reasonable way to overcome such a difhculty is to replace the white-noise fields 
with more general ones: this will be the subject of papers in preparation. Such 
attempts, however, require a preliminary and detailed investigation aimed to 
clarify that the consideration of noises which are non white in time docs not 
lead to inconsistencies and preserves the nice features of the models based on 
noises which are white. The analysis we are going to perform is also of interest, 
per se, for the study of stochastic dynamical equations in general, a subjet which 
has received a lot of attention in recent times. 

Sections 3 to 6 of this paper are devoted to a general discussion of non- 
white random differential equations in Hilbert space. We will show that the 
most important features of white-noise stochastic equations, in particular the 
desired reduction effects, hold also in the more general case, thus proving that 
the white-noise character of the stochastic processes is not an essential element 
of the dynamical reduction program. Dynamical reduction models based on 
general Gaussian noises have been studied by Pearle |^ and Diosi et al [|o[ ^ 
along different lines. In the final section of the paper, we apply the results and 
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the formalism of the previous sections to a specific model of dynamical reduction 
and we prove that it leads precisely to the localization of macroscopic objects 
in space. 



2 Review of CSL 

The CSL (Continuous Spontaneous Localizations) version of Dynamical Reduc- 
tion Models 1^, ^ is based on a stochastic Schrodinger equation which, in the 
Stratonovich language, takes the form: 



dm)) 

dt 



Ho + Y.A,w,{t) - 7^]^' 



\m)- (2.1) 



Here, Hq is the free Hamiltonian of the system; {Ai} is a set of commuting self- 
adjoint operators (representing the preferred basis) whose common eigenmani- 
folds are the linear manifolds into which the statevectors of individual physical 
systems are driven; Wi{t) are c~number independent stochastic processes with 
a gaussian distribution which is white in time: 

{{w,{t))) = 0, {{w,{h)wj{t2))) =-f5,j S{h - h), (2.2) 

the symbol ((.)) denoting the stochastic average associated to the process 



Equation (2.1) describes a diffusion process in Hilbert space; it is a linear equa- 
tion like the Schrodinger equation, but it does not preserve the norm of \'ip{t)) 
since the evolution is not unitary, due to the presence of the last two terms on 
the right hand side. The solution {ipit)) cannot therefore be endowed with a 
direct physical meaning. 

To overcome this difficulty, and at the same time to ensure that the reduc- 
tion mechanism reproduces the quantum mechanical probabilities, the following 
strategy has been adopted P, p[ . The physical vectors are the normalized solu- 



tions of equation (2.1) 



and it is assumed that any particular realization of the stochastic processes 
Wi{t), yielding the state \ipphys(t))^ has a probability of occurrence Pcoak[w{t)] 
equal to: 

PcMt)] = P^,.^,[wit)]mt))f, (2.4) 
where Pji^^[w{t)] is the original probability distribution of the gaussian white 
noises given by (|2.2|). 
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Of course, since ( |2.4| ) defines a probability distribution, it must sum to 1: 
V[wit)] Pc^^4w{t)] = 1 (2.5) 



= {((mm)))) 

(we remember that ((.)) refers to the average with respect to the original prob- 
ability distribution P„j^^[w{t)]). Equation ( pTs] ) imposes that the stochastic av- 
erage of the square norm of the vector {ipit)) be conserved; as it can be easily 



verified, equation (2.1) guarantees that this is the case. 



It is possible to write Q the norm-preserving equation for the physical vector 



IV'phy=(0) which is equivalent to eq. ( ^.iD with the prescription (2.4): such an 
equation is nonlinear and more difficult to handle. We will not consider it here. 
In reference Q it has been shown that, if one ignores the Hamiltonian term 



Hq, equation (2.1) together with the cooking prescription (2.4) — or, alterna- 
tively, the corresponding norm-preserving equation — drives the statevector 
of any individual physical system into one of the common eigenmanifolds of 



the operators At. Here we will simply show that equation (2J) implies the 
diagonalization of the density matrix with respect to the basis of the common 
eigenmanifolds of the operators^] Ai. 

The statistical operator is the average value, with respect to the the "cooked" 
(i.e. the physical) probability distribution Pcaak[w{t)], of the projections opera- 
tors onto the one dimensional linear manifolds spanned by the physical vectors 
IV'phy=(0) : 



pit) = y I?[u;(i)]|V'phy=(i)>(V'phy4i)l^coo.[u'(t)] 

= ((m))(m\))- (2.6) 

We note that, from the mathematical point of view, p{t) corresponds also to 
the ensemble of operators \ij{t))(ijj{t)\ {\ijj{t)) being non normalized), averaged 
with the raw probability distribution Pji^^[w{t)]. Thanks to this property the 
dynamical evolution equation for p(t) can be easily derived: 

^ = -l.[Ho,p{t)] - lY.^A.,,[A,pm. (2.7) 



^Of course, the diagonalization of the density matrix is only a nenessarv not, a sufficient 
condition for the localizations to occur, as it has been shown in reference [|l3|. Anyway, as 
already remarked, reference P| contains a proof that equation does imply the reduction 

of the statevector into the desired eigenmanifolds. 
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To show the effect of the reducing terms, let us suppose, for simplicity, that 
the common eigenmanifolds Ma of the operators Ai (which we assume to have 
a purely discrete spectrum) are one-dimensional and we call \a) the vector 
spanning Mq,: 

A,\a) = a,a\a). (2.8) 



For the moment, let us ignore the Hamiltonian term Hq. Then, equation (2.7) 
implies the following equation for the matrix elements {a.\p{t)\(3): 

d{a\p{t)\[3) ^ ( n2 / I (.\\R\ (on\ 
= ~ l^2^{aia - a.ip) {a\p{t)\l3). (2.9) 

i 

Equation ( |2.9| ) shows that the off-diagonal elements of the density matrix, cor- 
responding to the interference terms arising from the superpositions of different 
eigenstates of Ai, are exponentially damped. The diagonal elements, on the 
other hand, do not change with time. 

Up to now we have described the general formal structure of CSL. To give 
a physical content to the model, one must choose the "preferred basis" , i.e. the 
operators Ai which define the manifolds onto which the wavefunction is reduced. 
Obviously, our aim is to induce the spatial localization of macroscopic objects. 
To this purpose one can make the choice 

A, ^ AA(x) = (^)"^ /'d3ye"2^''"^\t(y^,)„(y^,)^ (2.10) 

where a^(y, s) and a(y, s) are the creation and annihilation operators for a con- 
stituent with spin component s, at point^ y. The parameter \j \/oi is a measure 
the localization accuracy of the reducing mechanism: for physical reasons |l[] its 
value has been chosen to be I/a/o — 10~^cm. The value of the other parameter 
of the theory, 7, which measures the strength of the correlation function of the 
white noises, is related to the parameter A ~ 10"^^ sec~^ of QMSL jl], which, 
within such a discrete model, specifies the frequency of the random reduction 
processes, according to 7 = A(47r/a)'^/^. Accordingly, the stochastic processes 
Wi{€) are replaced by a gaussian stochastic field ■u;(x, i), whose first two moments 
are: 

((«;(x,t))) =0, ((u;(x,ii)u;(y,<2))) =7 5(x-y)5(ti-t2). (2.11) 



■^If more than one type of particle is involved, an extra sum over the different kinds of 
particles must appear in the definition of A^(x). 
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The modified Schrodinger equation (2.1) becomes then: 



d\m) 



dt 




(2.12) 



and the corresponding equation for the statistical operator is: 




= -^[Ho,pit)] - I 1 d'x [AA(x),[AA(x),p(t)]] 



(2.13) 



Equations ( ^.12 ) and (2.13) exhibit two basic features: 



• At the microscopic level — i.e. when only few constituents are involved — 
the new terms do not alter in any appreciable way the pure Schrodinger 
evolution: all quantum properties of micro-systems are left essentially 
unchanged. This is due to the fact that the value of A is extremely small. 

• At the macroscopic level, on the other hand, the new terms induce in a 
very short time — much shorter than the perception time of a conscious 
observer — the suppression of the superposition of different macroscopic 
states and the reduction to one of them. Accordingly, macroscopic objects 
are always localized in space, and their classical properties are restored. 

This is the way in which dynamical reduction models are able to account for 
the behaviour of both microscopic quantum and macroscopic classical systems. 

3 Dynamical reduction models with general Ga- 
ussian noises 

In this section we begin the analysis of dynamical reduction models in which the 
reduction mechanism is controlled by general Gaussian noises. The first task 
is to derive a modified Schrodinger equation generalizing equation (^), and 
preserving the average value of the square norm of vectors, so that the cooking 
prescription can be applied to it. 

Let us then consider the following equation: 



where, as before, Hq is the Hamiltonian of the system, {Ai} is a set of commuting 
self-adjoint operators, and Wi{t) are c-number gaussian stochastic processes 




(3.1) 
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whose first two moments arcj^ 



(3.2) 



As in standard CSL, tlie evolution described by equation (3.1) is not unitary 
and it does not preserve the norm of the statevector; we then follow the same 
prescription outlined in section 2. We consider as physical vectors the normalized 
ones: 



II^W>II' 

and we assume that any particular realization of the stochastic processes 'Wi{t) 
has a probability of occurrence Pcook[w(t)] equal to: 



(3.4) 



where PB^^^[w{t)] is now the gaussian probability distribution defined by (3.2). 
The above assumptions guarantee, as we will show in section 5, that the reduc- 
tion probabilities reproduce standard quantum mechanical probabilities. 

As in section 2, we have to impose that equation (3.4) correctly defines a 
probability distribution, i.e. that it sums to 1. From equation (2.5) we see that 
this is equivalent to requiring that the time derivative of {{{'ip{t)\^l;{t)))) is zero. 
Let us evaluate it: 



dt 



mtmrn) = 




The two terms involving the Hamiltonian Hq cancel out (in fact they describe 
the unitary part of the evolution); the noises Wi{t), being c-numbers, can be 
taken out of the scalar product, so that: 
d 



dt 



mmm) = 2Y^{mt)\AM{t))^^{t))). 



(3.5) 



The right hand side of (3.5) can be rewritten with the help of the Furutsu 
Novikov formula : 

'5F[w{t)] 



{{F[w{t)]wm 



7E 

i 



Aj(i,s) 



dwj{s) 



ds 



(3.6) 



''There is no loss of generality in considering gaussian processes with zero mean. In fact, 
if ((uii(t))) = mi{t) 7^ 0, we can always define new proc esses Zi{t) = Wi{t) — mi{t), which have 
zero mean, and rewrite the Schrodinger equation (p.ll) in terms of the processes Zi{t). 
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(for simplicity, throughout this section we take io = as the initial time). 
F[w{t)] is any functional of the stochastic fields 'Wi{t); in the present case case, 

F[w{t)] = {m\Am))- 

The formal solution of equation (|3.lD is: 



Wi{s) ds 

\m)- 



(3.7) 



Note that, since \ip{t)) depends on the stochastic processes Wi{s) only within 
the time-interval [0,t], the functional derivative of \tp{t)) with respect to Wj{s) 
is zero if s ^ [0, <]. We then have: 



dt 



mtmm) 



27 E / ^»^(*'^) 



27 



1,3 



i,3 



mt)\ 



A,\^{t)) )) ds 



D,,it,s){{m)\A. 



sm)) 



Swj{s) 



(3.8) 
ds ^ 0. 



Since the time derivative of the average value of the square norm of the statevec- 
tor is not zero, we have to add an extra term to equation ( |3.lD , as expected and 
as it happens also in the case of white noise. Relation ( |3.§| ) tells us which kind of 
term must be added. The conclusion follows: with reference to our procedure, 
the request that Pcoak[w(t)] correctly defines a probability distribution, i.e. that 
the average value of the square norm of the statevector \ip{t)) is conserved, leads 
to the stochastic Schrodinger equation: 



d\m) 

dt 



Ho + ^A,w,(t)-2-i^A, [ dsD. 
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Swj{s) 



\m)- 



(3.9) 

This is the main result of this section. Note that an equation like (3.£) has been 



derived in references |l3| by following a different line of thought. 
Some comments are appropriate: 



• Equation ( |3.9| ) no longer describes a Markovian evolution for the state- 
vector unless the correlation functions Dij(t, s) are Dirac-J's in the time 
variable — i.e. the stochastic processes Wi{t) are white in time. As a 
consequence, the corresponding equation for the statistical operator is 
not of the quantum-dynamical-semigroup type |14|, n5| , contrary to what 



happen for the case of CSL (see equation (2.7)) 
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• In general, the explicit form of the functional derivatives of \^(t)) with re- 
spect to the noise Wi (t) cannot be evaluated exactly, except for few special 
cases, two of which will be considered in the next section. Therefore, in 
the general case it is difficult to analyze the time evolution of the statevec- 
tor and the statistical properties of the ensemble of states generated by 
the stochastic processes. In particular, one cannot write a closed equation 
for the evolution of the statistical operator. 



4 Two special cases 



In order to understand the kind of difficulties one encounters when working with 
non-white stochastic processes, and in particular the reasons for which the func- 
tional derivative of the statevector {ipit)) in general cannot be computed exactly, 
let us reconsider equation (^^, writing explicitly its perturbative expansion: 



Te 



oo 

= E 



Wi{s) ds 



i 


n 










f dh-- 


■L 






nl Jo 





f dt„T{H{ti)...H{t„)}, (4.1) 
Jo 



where we have defined the operator: 

H{t) = Ho + ih'^AiW^it). 



(4.2) 



The functional derivative of |'0(i)) with respect to Wj{s) can be obtained 
deriving term by term the series^ ( [4.l| ). The derivative of the term n = is 
zero; the derivative of the term n = 1 is: 



6wj{s) 

The next (n — 2) term is: 



if* 1 if* 

- / dtiH{ti) = -- / dti [inS{s-ti)Aj] = Aj 
"Jo J ^0 



i 













(4.3) 



(4.4) 



The functional derivative of the time-ordered product T{H{ti) H{t2)} = 0{ti 
t2)H{h)H{t2) + e{t2 - h)H{t2)H{ti) is: 
<5 



dwj{s) 



T{Hih)H{h)} 



= ihe{h-t2)[5{h- s)AjH{t2) + S{t2- s)H{h)Aj] 

+ ihe{t2-h)[5{t2-s)AjH{ti) + 5{h-s)H{t2)A,]. (4.5) 



We assume that the initial state |</'(0)) does not depend on the stochastic processes Wi{t). 
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We note that the first and third terms at the right hand side of (4.5) differ only 
for the exchange of the dummy variables ti ^2 ; the same is true for the second 
and the fourth term. The derivative of the 71 = 2 term (i.e. of Eq. (4.4)) is 
then: 



A, 



dti H(ti) 



A. 



(4.6) 



Equation (4.6) does not have a simple form, contrary to ( |4.3| ), and derivatives of 
higher terms are more and more complicated, due to the fact that the operators 
Aj in general do not commute with the Hamiltonian Hq. In fact, would they 
commute, equation (4.6) would simplify to: 

A, 



dti H{ti 



(4.7) 



i.e. the derivative of the second term would give Aj times the first term. More- 
over, if [Aj, Hf)] ~ 0, the functional derivative of the term n + I gives Aj times 
the n-th term: 



Swj{s) 



\m) 



A, \m), 



(4.8) 



as we are going to prove. In fact, the hypothesis that the operators Ai commute 
with the Hamiltonian Hq is equivalent to the (more elegant) requirement that 
the operators H{t) defined in (|4.2|) commute at different times. In this case, the 



time-ordered product in the exponential series (4.1) can be omitted, and the 
functional derivative of the n-th term is: 



5wj{s) 



i 


n 










- f dti-- 


■I 






"! Jo 





A, 



1 " 

1=1 
1 

(n-1)! 

-1 

(n 







f dti-- 









dt,, {H{ti) . 



H{ti 



■ H{tn)} = 
SHiU 



6 Wj{s) 



■ H{t, 







f dti-- 









SH{ti 



H(tn) 



5 Wj (s) 

i— / dti--- j dtn-i {H{ti)...H{tn-i)}. (4.9) 
~ ^)- Jo Ja 



This completes the proof. Note also that, when s = t, an extra factor 1/2 
appears in (4.S), because in this case the Dirac delta function arising from the 
functional derivative of H{t) is centered in one of the two extreme points of the 
interval of integration. 

Recently, S. Adler and P. Horwitz ||l^ (see also ||l^) have proposed a white- 
noise model of dynamical reductions in which the operators Ai are taken to be 
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functions of the Hamiltonian Hq; this implies that the stochastic terms of equa- 
tion drive the statevector into the energy eigenmanifolds of the physical 



system. Making such a choice in the non- white equation (3.9), the operators 
H{t) at different times commute among themselves, the functional derivatives 



of the statevector \ip{t)) can be computed, and equation (3.9) becomes: 



d\m) 

dt 



E 



Aiw^ (t) - 27 ^ AiAj / D,j {t,s)ds 



1,3 



\m), (4.10) 



with Ai = Aii^Ho). Equation ( 4.1C ) is exact and, correspondingly, one can easily 
derive a closed equation for the time evolution of the statistical operator. All the 
statistical properties concerning the physical system can be evaluated exactly. 

We conclude the section showing that the functional derivatives of \ip{t)) 
can be explicitly evaluated also in the case of general white noise stochastic 
processes, without having to require that Hq commutes with Ai. Moreover, we 



will prove that in this case equation (3.9) reduces to (2.1), as expected. 

Under the assumption of white-noise stochastic processes {Dij{ti,t2) 
Sij 5{ti — ^2)), the Furutsu-Novikov relation 



{{F[w{t)]w,{t))) = 7 



mwjt)] \\ 

Sw,{t) 11 



(4.11) 



leads to the following expression for the time derivative of the average value of 
the square norm of the statevector |'0(t)) satisfying equation (|3.lD: 



dt 



immm) = 27^ 



_ Swi{t) 



A^\m) 

Swi{t) 



(4.12) 



We now have to evaluate the functional derivatives of the statevector, taking 
into account that the noises Wi{t) (appearing in the derivatives) are taken at 
time t. 

The derivative of the term n = 1 is equal to {1/2) Aj (see equation (^T^)), the 
factor (1/2) deriving from the Dirac delta function d{t — ti) which is integrated 
between and t. For the derivative of the n = 2 term, let us look at expression 
(^^). If we take s — t, the second term goes to zero, while the first one gives^: 

ct 



2 ■' 



dti H{ti 



(4.13) 



^The factor (1/2) appears for the same reason as before. 
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In general, the functional derivative of any terms of the exponential series (4.1) 
gives (1/2) times the previous term, so that: 

S 



1 



This means that the square-norm-preserving Schrodinger equation is: 

d\m) 



dt 



(4.14) 



(4.15) 



which coincides with the original CSL equation (2.1). An alternative and quicker 



way to derive the white-noise limit is to replace Dij(t,s) with 6ijS{t — s) in 



equation (3.9) and to show that (4.14) is a consistent solution 



5 The reduction mechanism 



Here, we will analyze under which conditions the new terms in the modified 



Schrodinger equation (3.9) induce, for large times, the reduction of the stat- 
evector to one of the common eigenstates of the commuting operators Ai. 

For this purpose, let us disregard the Hamiltonian Hq; under this assumption 
the operators H{t) commute at different times and (as discussed in the previous 
section) the functional derivatives of the statevector \ip{t)) give the operators 
Ai times |'0(^))- Equation ( |3.9| ) becomes therj^: 



dm)) 

dt 



A, (t) - 27 ^ A,Aj I D,j {t,s)ds 



\m)- 



(5.1) 



The equation for the statistical operator can now be easily derived; using the 



definition ( |2.6|) , we get: 

dpjt) 
dt 



-jY.[A,,[Aj,p{t)]] f D,,{t,s)ds, 

•Itn 



(5.2) 



which is a consistent generalization of the CSL equation (2.7) when the Hamil- 
tonian iJg is omitted: in fact, if the stochastic processes 'Wi{t) are independent 
and white (Dy(<i,t2) = 5ij S{ti — ^2)), then (5_^) reduces exactly to (2.7). 



In order to test the reduction properties, we will show first of all how the 
reduction mechanism works for the statistical operator (see footnote 1). As in 
section 2, let us suppose that the common eigenmanifolds of the operators A^, 

^Here and in what follows, we consider a generic initial time to. 
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which we assume to have a purely discrete spectrum, are one-dimensional; let 
I a) be the vector spanning the a-eigenmanifold. The equation for the matrix 
elements {a\p{t)\P) is: 

d{a\p{t)\l3) ^ -^^{a^a- aip){aja- ajf}) f D.i-j{t, s) ds {a\p(t)\l3) . (5.3) 

Making use of the symmetry property of the correlation functions: 

D.jitiM) = Dj,{t2,ti), (5.4) 



we can write the solution of equation (5.3) in the following form (see also ]19|): 



- -^^{uia - aip){aja - ajp) / dti / (it2 (ti, ^2) 
{a\pitm=e («|p(to)|/3). 

(5.5) 



From equation (5.5), we sees that if \a) — \(3), the exponent is zero: as in CSL, 
the diagonal elements of the density matrix do not change in time. If, on other 
other hand |a) ^ \(3), the evolution of the matrix element depends on the time 
behavior the correlation functions Dij(ti,t2). 

If we want the off-diagonal elements to be damped at large times, two con- 
ditions must be satisfied. The first one is that the exponent in (5^) must 
be negative: this is always true, since the correlation function of a Gaussian 
process is positive definite. 

The second condition is that the double integral of the correlation function 
must diverge for large times: 

/ dti / dt2 Dij{ti,t2) — > +00 for t ^ +CX), (5.6) 
J to J to 

so that the off-diagonal elements of the density matrix go to zero. This condition 
is not a priori satisfied by a generic Gaussian stochastic field. At any rate, 
physical reasonable stochastic fields always satisfy it: here we present just a 
couple of meaningful examples. 

Suppose the stochastic fields Wi{t) are equal and independent, with a (nor- 
malized) Gaussian correlation function: 

D^jihM) = 2r2 . (5.7) 



Let us also take to = —00. Equation ( |5.3| ) then becomes: 
d{a\p{t)\P) 7 



dt 



Y,{a,o.-a,pf{a\p{tm, (5.8) 
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which is independent from the correlation time r, and moreover it corresponds 
exactly to the CSL equation (2.£). The fact that we have taken to — —oo means 
that the correspondence between equation (|5.3|) — with a correlation function 



like (5.7) — and equation (2.9) is exact only in the limit of large times. Note 
also that if we take the limit r ^ 0, the gaussian process becomes a white 
noise process with a Dirac-J correlation function and we recover, again, the 
CSL theory 

As a second example, suppose the correlation function is: 



Dij{ti, t2) 

Equation (^) becomes: 

d{a\p{tm _ 



\ti-h\ 



(5.9) 



dt 



1 - e 



^(a»a-a,^)2(a|p(t)|/3). (5.10) 



As before, the off-diagonal elements are exponentially damped and, in the limit 
t +00 we recover the behavior of CSL. Note that the effect of a non-white 
correlation function is that of decreasing the reduction rate of the localization 
mechanism. 

We now analyze how the reduction mechanism works at the wavefunction 



level, proving in this way that equation (5.1) leads to the reduction of the 
statevector into one of the common eigenmanifolds of the operators Ai. As 
in reference we consider a simplified dynamics in which only one operator 
A appears in equation This operator is coupled to a single stochastic 

process w(t), whose correlation function is D[ti, t2). Finally, we assume that at 
the initial time Iq the statevector is: 



(5.11) 



where Pa and Pp are projection operators onto the eigenmanifolds of A cor- 
responding to two different eigenvalues a and /3, respectively. The solution of 
equation (5.1) is: 



= e"^W-"V(i)p„|^(to)) + e/3^W-/3V(i)p„|^(io)), (5.12) 

where 



x{t) 



w{s) ds, 



fit) 



I dsi f ds2D{si,S2). (5.13) 

J to Jto 
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Note that 7/(t) = {{x'^{t))), i.e. such a quantity is the variance of the stochastic 
process x{t). 

Since the "raw" probabihty distribution of the process x{t) is: 



P..^[x{t)] = -==e 27/W , (5.14) 
\/27r7/(0 



taking into account the cooking prescription (3^) we obtain: 

^ 'x{t) - 2/37/(t)]2 



(5.15) 



Equation (5.15) imphes that, if f{t) +00 when t —>■ +00, the stochastic 
process x(t) wiU take either a value close to 2a"ff{t) — within an interval of 
width \/"ff{t) — or a value close to 2(3^ fit), within the same interva][]. Of 
course, the requirement that f{t) +00 as time increases is exactly the same 



as requirement (5.6) which guarantees the damping of the off-diagonal elements 
of the density matrix. 

Suppose now that the actual realization of the stochastic process x{t) occurs 
around 2ajf{t); the corresponding probability is ||Pct|7/'(0))|p. We then have: 



wmmn' ^ -2Ma-0?f(t) \\mm)\\ 
WPcimw ~ \\Pa\m)\\ 



2 



sls t ^ 00, (5.16) 



which means that the statevector \ip{t)) is driven into the eigenmanifold of the 
operator A corresponding to the eigenvalue a. By the same reasoning, it is 
immediate to see that, with a probability equal to ||P/3|i/'(0))|p, the statevector 
is driven into the eigenmanifold associated to the eigenvalue /?. We have thus 
proved that the statevector \ip{t)) undergoes a random spontaneous localiza- 
tion into one of the two eigenmanifolds of the operator A, with a probability 
which coincides with the one assigned by standard Quantum Mechanics to the 
outcomes of an experiment aimed to measure the observable A. 

^As noted in even though the interval if if) tends to infinity as time increases, the 
ratio ^ 7/(i)/2(o — f3)-ff{t) goes to zero. 
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6 The average value of observables 



When one disregards the Hamiltonian term Hq, it is not difficult to see how the 
stochastic terms affect the average value of physical quantities. 

The mean value of an operator O (for simplicity, we consider an observable 
which does not depend explicitly on time) is defined as the expectation value 
{(j){t)\0\(j){t)) , averaged over all possible realizations of the stochastic noises: 



(o) = J v[w{t)]m)\o\m)Pc...[wit)] = 
= mt)\o\m)))- 

Its time derivative can be calculated following almost the same steps which, in 



the previous section, have led to equation (5.2) for statistical operator; the final 
equation is: 



to 



to be compared with the corresponding CSL-white noise equation: 

^ = -|^(((^(t)|[A.,[A„0]]|V.(t)»). (6.2) 

i 

The analysis of the previous section should have made clear how (6.1) differs 



from (6.2), so we will not repeat it here. 



7 Connection with CSL 

We now apply the formalism introduced in the previous sections to derive an 
equation with the property of localizing macroscopic systems in space, like in 
CSL; in other words, we specify the choice of the "preferred basis" {Ai} in such 
a way to have a physically meaningful theory for our purposes. 

The most natural choice for the operators Ai is the number density operator 
for a system of identical particles: 

a\-x, s) a{-x, s) . (7.1) 

s 

Correspondingly, the noises Wi{t) are replaced by a stochastic field w(x, t), whose 
correlation function is -D(x, ii; y, 12). 

In reference Q], the transformation and invariance properties of dynami- 
cal reduction models have been discussed in detail. In particular, it has been 
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proved that, in order for the physical properties of the model to be invariant 
under Galilean transformations (we speak of stochastic Galilean invariance), the 
correlation function D(x, ii; y, ^2) itself must be invariant under the considered 
group of transformations, i.e. 

Z?(x,ti;y,i2) = D(|x-y|,ti -i2); (7.2) 



the easiest way to construct a function like (7^) is to take the product of two 



functions of the space and time variables, respectively: 

Z?(x,ti;y,i2) = g{\x-y\)h{ti-t2). (7.3) 
As regards g(|x — y|), a reasonable choice is a gaussian function, like in CSL: 

5(|x-y|) = 7(£)' n^^"^^ , (7.4) 

with 1/y/a ~ 10~^ cm. 

It is natural to choose a gaussian function also for h(ti — t2): 

/ /? \ 5 -^(ii - t2)2 

h{t,-t,) = i^f-j e 4^' " . (7.5) 

With the above choice, we have introduced a new parameter (/?); this can be 
considered as a drawback of the model. However, we note that it always is 
possible to define (3 in terms of a, 7 and fundamental constants of nature, so 
that no new arbitrary parameter is introduced into the model. As an example, 
we can choose /? = c^a ~ 10'^° sec~^, where c is the speed of light. This choice 
is particularly appropriate in the light of a possible relativistic generalization 
of the theory, which we will discuss in a future paper. Moreover, such a choice 
corresponds to an extremely small correlation time, so that for ordinary systems 
(moving slower than the speed of light) the behavior of the model is similar to 
the one deriving from the white-noise CSL. 

The modified equation ( |3.9| ) for the statevector evolution becomes now: 

d\m) 

dt 



-^-ffo + y d^a:7V(x)w(x,i) - (7.6) 
-27y d'xd'yN{^)g{\^~y\) dsh{t-s) ^^ |V^(i)). 
If we ignore the free Hamiltonian Hq, i.e. we confine our considerations to the 
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reduction mechanism^, equation (7.6) becomes 

dim) ^ 

dt 
with: 



d^a;iV(x)u.(x,<) - ^{t) / d'xd'y N{^)g{\^ - y\)N{y) 



m)) 

(7.7) 

7(t) = 27 / dsh{t^s). (7.8) 

Jto 

The corresponding equation for the statistical operator is: 

|p(t) = 1 d3^rf3y[iV(x),[iV(y),p(t)]].9(|x-y|). (7.9) 



Equation (7.6) can be rewritten in a form closer to equation (2.12), which 
will be useful for the subsequent discussion. Let us define a new Gaussian 
stochastic process ui(x, i), which is connected to w(x, t) by the relation: 



gy Id^.e 2 



(x-y)^ 



w{y,t). 



w(x, t) = 

Iy(x, t) has zero mean, and correlation function 

((W(x,ii)w(y,t2))) = 7'5('Hx-y)/»(ti-t2). 



(7.10) 



(7.11) 



Using the following relation: 



Sw{x, s) 



(5w(x, s) Sw{y, s) 



a \2 
2^ 



d-'ye 2 



Sw{y,s) 



it can be easily seen that (7.6) is equivalent to the equation: 

d\m) 

dt 



-^iJo + / d^xAf(x)w(x,t) 



(7.13) 



27 / d^xM{x.) / dsh{t-s) 



to 



Sw{x, s) 



\m), 



with Af{x) defined by ( |2lo| ). 

*For the physically interesting cases, e.g. for the dynamical evolution of macrosystems, 
such an assumption is justified by the fact that the effect of the reduction is much faster that 
the tipical times in which the Hamiltonian can induce appreciable dynamical changes of the 
statevector. 
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7.1 Dynamics for macroscopic rigid bodies 

As for white-noise CSL |3|, it is not difRcult to discuss the physical imphcations 



of equation ( |7.q ) — or equation (7.13) — for the case of a macroscopic rigid body, 
i.e. a body such that the wavefunctions of its constituents can be considered 
very well localized with respect to the localization length 1 / 

To be precise, in analogy with the procedure followed in , let us consider 
a system of IS! identical particles of coordinates qi; let 



1 ^ 



be the center of mass coordinate, and let us write 

q,: = Q + q^, (7.15) 

where the coordinates qj are functions of 3A^ — 3 independent internal variable^, 
which we call r. Let us consider the wavefunction 



|V'(g,s)) = |0(Q))|^(r,s)) |^(r,s)) 



A 
S 



\A{r,s)), (7.16) 



where q — {qi} and s — {si} are the sets of the space and spin coordinates 
of the N particles, respectively, while "A" and "S" mean symmetrization or 
antisymmetrization with respect to the interchange of the variables (qi,Si), 
respectively. 

In reference |^ it has been proved that if the wavefunction of the internal 
degrees of freedom is very well peaked with respect to the characteristic length 
l/-y/a, then, to an extremely high degree of accuracy, 

AA(x)|0(Q))|^(r,s)) = FiQ~^)\m))Mr,s)), (7.17) 

with 

i=l 

where rp describes the set of the average equilibrium positions of the particles 
of the rigid body. Equation (7.18) means that the operators A/'(x) act only on 
the center of mass wavefunction |0(Q)). 

As a consequence, if the Hamiltonian Hq can be written as 

= Hq + Hr, (7.19) 

^See ref. |^ for further details on the degrees of freedom of the system. 
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and if |(^(Q)) and Ifir, s)) satisfy the equations 

d\m,t)) 



dt 



h 

27 



Q 



d^x7V(x)w(x, t) 



(7.20) 



d^xj\f{x) / dsh{t-s) — 



(5w(x, s) 



d\ip{r,s,t)) 
dt 



-Hr 



|(/?(r,s,t)), 



(7.21) 



then \ip{q,s,t)) satisfies equation ( 7.13| ) or, equivalently, equation (7.6). 

Equations ( [7.2C| ) and ( 7.21 ) imply that the center of mass and internal mo- 
tion decouple, and that the stochastic terms affect only the center of mass and 
not the internal structure, as it happens for CSL. 

Following the same arguments of reference H), it can also be proven that 
the localization rate of the center of mass wavefunction grows linearly with 
the number of particles of the rigid body. Such a localization rate can be 
easily computed by studying the time evolution of the off-diagonal elements 
(Q'I/OqIQ") of the reduced statistical operator describing the center of mass 
motion of the system. As we did in section 5, we disregard the Hamiltonian 
iJq, in accordance with the fact that, for |Q' — Q"| > the reduction rate 

turns out to be much faster than the typical times required in order that the 
standard Schrodinger evolution induces appreciable changes of |0(Q, t)). Under 
this assumption, equation ( 7.2C ) becomes: 

d\m,t)) 



dt 



d^a;i^(Q-x)w(x,t)-7(t) /d^a;F^(Q-x) 



l0(Q,t)), 

(7.22) 



and the corresponding equation for the matrix elements (Q'IpqIQ") is: 



d(Q'|pQ(t)|Q^' 
dt 



-r(Q',Q",i) (Q'lpQ(t)IQ" 



(7.23) 



with 



r(Q',Q",t) = 7(t) j d 

F(Q' - x)i^(Q" - X 



^^'(Q'-x) 



V(Q"-x) 



(7.24) 



This is the same term as the one appearing in CSL, with ^{t) replacing 7; this 
proves that also in our model the reduction frequency of the center of mass of 
the system grows linearly with the number of its constituents. Moreover, taking 
a large value for (3, as it has been suggested previously, ^{t) — > 7 in very short 
times, so that the reducing dynamics is practically the same as the one of CSL. 
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